Superfluid density and specific heat within self-consistent scheme for two-band 

superconductor 
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The two gaps in a two-band clean s-wave superconductor are evaluated self-consistently within 
the quasiclassical Eilenberger weak-coupling formalism with two in-band and one inter-band pairing 
potentials. Superfluid density, free energy and specific heat are given in the form amenable for fitting 
the experimental data. Well-known two-band MgB2 and VsSi superconductors are used to test the 
developed approach. The pairing potentials obtained from the fit of the superfluid density data in 
MgB2 crystal were used to calculate temperature-dependent specific heat, C(T). The calculated 
C(T) compares very well with the experimental data. Advantages and validity of this, what we 
call the "7- model", are discussed and compared with the commonly used empirical (and not self- 
consistent) "a-model". Correlation between the sign of the inter-band coupling and the signs of 
the two order parameters is discussed. Suppression of the critical temperature by the inter-band 
scattering is evaluated and shown to be severe for the inter-band repulsion as compared to the 
attraction. The data on a strong T c suppression in MgB2 crystals by impurities suggest that the 
order parameters on two effective bands of this material may have opposite signs, i.e., may have the 
s± structure similar to the current proposals in iron-based pnictide superconductors. 

PACS numbers: 74.25.Nf,74.20.Rp,74.20.Mn 



I. INTRODUCTION 

Nearly all superconductors discovered recently and 
some well-studied compounds (e.g., Nl^Se and VaSi) are 
multiband materials with complex Fermi surfaces and 
anisotropic order parameters. Measuring temperature 
dependences of the London penetration depth A(T), often 
converted to the superfluid density p(T) = A(0) 2 /A(T) 2 , 
and of the electronic specific heat C(T), are among pri- 
mary tests directly linked to pairing mechanisms of new 
superconductors. Still, the methods employed to inter- 
pret the data are often empiric with simplicity as a main 
justification. 

The most popular among practitioners "a-model" 
takes a shortcut by assigning the BCS temperature de- 
pendence to the two gaps Ai^ with which to fit data 
on the specific heat 1 and the superfluid density p — 
xpi + (1 — x)p2- 2,3 Here, p\ ,2 are evaluated with Ai^ = 
(ai,2/l-76)Asc , s(T) and x takes account the relative 
band contributions. Although the a-model had played 
an important and timely role in providing convincing ev- 
idence for the two-gap superconductivity in Mgl^, 1,2 it 
is intrinsically inconsistent in the most important task of 
this procedure, namely, in describing properly the tem- 
perature dependencies of p(T) and C(T). In fact, one 
cannot a priory assume temperature dependencies for the 
gaps in the presence of however weak interband coupling 
imposing the same T c for two bands. In unlikely situ- 
ation of zero interband coupling, two gaps would have 
single-gap BCS-like T— dependencies, but will have two 
different transition temperatures, see Fig. 1. 

The full-blown microscopic approach based on the 
Eliashberg theory, on the other hand, is quite involved 
and difficult for the data analysis. 4,5 Hence, the need for 
a relatively simple but justifiable, self-consistent and ef- 



fective scheme experimentalists could employ. The weak- 
coupling model is such a scheme. Over the years, the 
weak-coupling theory had proven to describe well multi- 
tude of superconducting phenomena. Similar to the weak 
coupling is the "renormalized BCS" model of Ref. 6 that 
incorporates the Eliashberg corrections in the effective 
coupling constants in a manner described below. We will 
call our approach a "weak-coupling two-band scheme" 
and clarify in the text below that applicability of the 
model for the analysis of the superfluid density and spe- 
cific heat data is broader than the traditional weak cou- 
pling. 

The s-wave weak-coupling multigap model has been 
proposed at the dawn of superconductivity theory by 
Moskalcnko 7 and Suhl, Matthias, and Walker 8 when nu- 
merical tools were still in infancy. In this work we basi- 
cally follow these seminal publications to develop a self- 
consistent procedure for the penetration depth data anal- 
ysis. Our scheme allows one to connect between two inde- 
pendent data sets, the superfluid density and the specific 
heat, thus providing a reliability check upon values of the 
coupling constants extracted from fitting the data. We 
also discuss the suppression of the critical temperature 
T c by non-magnetic impurities and suggest that the data 
on this suppression for MgB2 are consistent with a weak 
repulsive interband interaction that corresponds to op- 
posite signs of the order parameter on two bands, i.e., to 
± s structure of the order parameter. 

To test our formal scheme, the data for two known 
s-wave two-gap superconductors, MgB2 and VaSi, were 
used. The specific heat data were taken from the Ref. 9. 
The magnetic penetration depth, A(T), was measured 
by using a self-oscillating tunnel-diode resonator (TDR) 
with resonant frequency fo m 14 MHz. The measured 
quantity is the shift of this frequency, f(T) — fo = 



2 



—4:Ttx(T)Q , where \ 1S the total magnetic susceptibility 
in the Meissner state and Q ~ foV s /2V c (1 — iV) is a ge- 
ometric factor defined by volumes of the coil, V Cl and of 
the sample, V s , and by the demagnetization factor N. Q 
is measured directly by pulling the sample out of the coil 
at the lowest temperature. 3 For the susceptibility we use 
4ttx (T) = X(T)/w tanh[w/A(T)] - l. 10 



II. EILENBERGER TWO-BAND SCHEME 

Perhaps, the simplest formally weak-coupling ap- 
proach is based on the Eilenberger quasiclassical formula- 
tion of the superconductivity valid for general anisotropic 
order parameters and Fermi surfaces. 11 Eilenberger func- 
tions /, g for clean materials in equilibrium obey the sys- 
tem: 



= 2Ag/h-2ojf, 

a 2 = i-/ 2 , 

A(fe) = 2TrTN(0)^(y(k,k')f(k' 

w>0 



(1) 
(2) 

(3) 



Here, fc is the Fermi momentum; A is the gap function 
which may depend on the position k at the Fermi sur- 
face in cases other than the isotropic s-wave. Further, 
N(0) is the total density of states at the Fermi level 
per one spin; the Matsubara frequencies are defined by 
hu = nT(2n + 1) with an integer n, and lod is the Dcbye 
frequency; (...) stands for averages over the Fermi surface. 
As a weak coupling theory the Eilenberger scheme deals 
with the effective electron- electron coupling V responsible 
for superconductivity; properties of intermediate bosons 
(phonons or other possible mediators) enter via properly 
renormalized V. 

Consider a model material with the gap given by 



A(k) = Ai, 2 , k e Fi, 2 , 



(4) 



where Fi,F 2 are two sheets of the Fermi surface. The 
gaps are constant at each band. Denoting the densities 
of states on the two parts as -/Vi j2 , we have for a quantity 
X constant at each Fermi sheet: 

(X) = (XiNi + X 2 N 2 )/N(0) = ni Xi + n 2 X 2 , (5) 

where ni j2 = Ni, 2 /N(0); clearly, m + n 2 = 1. 

Equations (1) and (2) are easily solved. Within the 
two-band model we have: 

f* = ir> 9 » = ^> (% = Ai + nW, (6) 

Pv Pv 

where the band index v = 1,2. The self-consistency equa- 
tion (3) takes the form: 



a v = 2-kt ^2 n v x »»f» = X! ? v ; v a m ^2 ~o — > ( 7 ) 

(U=l,2 n u> ^ 



where X vtl = N(0)V(v, /z) are dimensionless effective in- 
teraction constants. 

A remark is here in order on applicability of Eq. (7) 
central for our approach. Starting with the general 
Eliashbcrg formalism, Nicol and Carbotte derived a 
"renormalized BCS" self-consistency equation: 6 



a, = (ash -/*:„)/„, 

^v 



(8) 



where Xf? is the coupling due to electron-phonon in- 
teraction, describes the Coulomb interaction, and 
Z v = 1 + XfJ^ is the strong coupling renormaliza- 
tion. Replacing (A^ — jZ v — » n„A„ M we obtain our 
Eq. (7). One should have this in mind while interpreting 
the constants A„ M which can be obtained from fitting the 
data to our "renormalized weak-coupling" model. 

Note that the notation commonly used in literature for 
Xvffl differs from ours: X„^ = n^X u ^. We find our nota- 
tion convenient since, being proportional to the coupling 
potential, our coupling matrix is symmetric: X vfl = A M „. 
It is worth stressing that for a given coupling matrix A M „ , 
relative densities of states n v , and the energy scale Hlod, 
Eq. (7) determines both T c and Ai >2 (T). 



A. Critical temperature T c 

As T — > T c , A12 — > 0, and (3 — > fko. The sum over to 
in Eq. (7) is readily evaluated: 

^ 2irT 2Huj d 2Huj d 

S = } — — =ln— =ln— — — , (9) 

^ huj T c ne-y 1.76 T c ' v ' 

7 = 0.577 is the Euler constant. This relation can also 
be written as 

1.76 T c = 2hw D e- s . (10) 
The system (7) is linear and homogeneous: 



Ai = 5(niAnAi + n 2 Ai 2 A 2 ) , 
A 2 = S(ni Ai 2 Ai + n 2 A 2 2A 2 ). 



(11) 



It has nontrivial solutions Ai j2 if its determinant is zero: 
S 2 n\n 2 r] - S(ni An + n 2 A 22 ) + 1 = 0, 



V — AnA 22 — Xf 2 . 
The roots of this equation are: 



(12) 



_ niAn + n 2 A 22 ± (uiXTi + n 2 X 22 ) 2 - 4nin 2 rj 
2n\n 2 r\ 

niAn + n 2 A 22 ± yj {n^XTi - n 2 X 22 ) 2 + ^n x n 2 X\ 2 



2n\n 2 r) 



(13) 



3 



Since T c <C hwr>, Eq. (10) shows that only positive S 
are admissible. The second form (13) shows that both 
roots are real; their product is S1S2 = \jn\nir\. Choos- 
ing a proper root, one should consider various possibil- 
ities. If 77 > (that implies both An and A22 are of 
the same sign), Eq. (13) shows that both roots are pos- 
itive, and one should choose the smallest (to have max- 
imum T c ). If 77 < (Af 2 > A11A22 that can happen 
(i) for a sufficiently strong interband coupling for both 
An and A22 positive or (ii) if one of An,A22 is repul- 
sive), one should take the square root with minus. It 
is of interest to note that even for An = A22 = 0, the 
interband coupling of either sign may lead to supercon- 
ductivity. In fact, S = 1/^/711712 I A12 1 for a dominant 
interband interaction | A12 1 3> l^iAn + U2A22|- However 
exotic, this possibility should not be ignored. This situ- 
ation has been considered time ago by Geilikman 12 who 
found that interband Coulomb repulsion could lead to 
superconductivity; recently this possibility has been con- 
sidered by Mazin and Schmalian in a discussion of su- 
perconductivity in the iron- pnictides. 13 If 7/ = 0, Eq. (12) 
yields S = 1/(tiiAh + 7i2A22)- Finally, if the interband 
coupling is exactly zero, a quite unlikely situation, the 
second form of S in Eq. (13) gives two roots 1/niAn and 
1/ti2A22- The smallest one gives T c , whereas the other 
corresponds to the temperature at which the small gap 
turns zero. This situation is depicted in Fig. 1. 

We conclude this incomplete list of possibilities by not- 
ing that within this model, interband coupling enters S 
only as \\ 2 , he., interband attraction in clean materials 
affects T c exactly as does the repulsion. This is no longer 
true in the presence of interband scattering, the question 
discussed below. Denoting the properly chosen root as 
5 = 1/A we have: 

1.76T c = 2/ia; £ >exp(-l/A). (14) 

One easily checks that for all A's equal this yields the 
standard BCS result. Among various possibilities we 
mention here the case 77 = A11A22 — A^ 2 = for which 

A = 711A11 + 712A22 = (A) . (15) 

This case corresponds to a popular model with factoriz- 
able coupling potential V(k,k') = V O,(k)fl(k') . 14 This 
potential is amenable for analytic work, but it reduces 
severely the richness of the two-band scheme. 

Since the determinant of the system (11) is zero, the 
two equations are equivalent and give near T c : 

A 2 A-TiiAn 

"T = \ ■ (16) 

Ai 712A12 

When the right-hand side is negative, A's are of oppo- 
site signs. Within the one-band BCS, the sign of A is 
a matter of convenience; in fact for one band the self- 
consistency equation determines only | A| . For two bands, 
Ai and A 2 may have opposite signs. If the A values are 
+Di and — D2, Eq. (7) shows that —D\ and +D2 is a 
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FIG. 1: (Color online) Calculated superfluid density and the 
gaps vs. reduced temperature (inset) for zero interband cou- 
pling, A12 = 0. In this calculation An = 0.5, A22 = 0.45, 
nl = n2 = 0.5, and 7 = 0.5, see text below. 
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FIG. 2: (Color online) Calculated superfluid density and the 
gaps vs. reduced temperature (inset) for A12 = 0.01, An = 
0.5, A22 = 0.45, nl — n2 = 0.5, and 7 = 0.5, see text below. 

solution too. One should be aware of this multiplicity of 
solutions when solving the system (7) numerically. The 
problem is even worse because Ai = A2 = is always a 
solution. 



B. Order parameter 

Turning to evaluation of A„(T), we note that the sum 
in Eq. (7) is logarithmically divergent. To deal with this 
difficulty, we employ Eilenberger's idea of replacing fkOr> 
with the measurable T c . These are related by Eq. (14) 
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which can be written as 



P = A 2 fc (0)/A 2 h (T): 



1 T ^2nT 



(17) 



A T c ^ hu 

UJ 

Now add and subtract the last sum from one in Eq. (7): 



= E n^Xy^n 



E 



/2ttT 2ttT 



E 



2ttT 



^ /2ttT 27tT\ 1 , T 



(18) 



The last sum over uj is fast-converging and one can re- 
place ud with oo. Numerically, the upper limit of sum- 
mation over n can be set as a few hundreds that suffices 
even for low temperatures. Introducing dimensionless 
quantities 



2-kT 



T r 2irt 



(19) 



with t = T/T c , we rewrite Eq. (18): 



^ = E "m^A ( \ + ln TF - A n) 

u=i,2 \a j y 



^ = E 



\ n + X / 2 ^ + (n + 1/2)2 



(20) 



For given coupling constants A„ M and densities of states 
n v , this system of two equations can be solved nu- 
merically for b v and therefore provide the gaps A„ = 
2wT5v(t). Two examples of these solutions with the sets 
of parameters differing only in A12 are given in insets to 
Figs. 1 and 2. We observe that even a small interband 
coupling changes drastically the behavior of the small 
gap. 



C. Superfluid density 

Having formulated the way to evaluate A(T), we 
turn to the London penetration depth given for general 
anisotropies of the Fermi surface and of A by (see e.g., 
Rcf. 15): 



167r 2 e 2 iV(0)r 



ik 



E 



P 3 



(21) 



where Vi is the Fermi velocity. We consider here only 
the case of currents in the ab plain of uniaxial or cu- 
bic materials with two separate Fermi surface sheets, for 
which a simple algebra gives for the superfluid density 



P 

Pv 

7 



7Pi + (1 - i)p2, 

00 

^]T[^ + (n+l/2) 2 ]- 3/2 



n=0 



(22) 



where v 2 , are averages over corresponding band of the 
in-plane Fermi velocities. The formal similarity of the 
first line here to the widely used a-model (with an x 
in place of out 7) prompts to name our scheme the "7- 
model" (we do not renormalize the gaps, so no analog 
of a exists in our scheme. We note, however, that our 
7 that determines partial contributions from each band 
is not just a partial density of states m of the a-model, 
instead it involves the band's Fermi velocities. 

We now apply the approach developed to fit the data 
for the superfluid density of MgB 2 crystals acquired by 
using the TDR technique described above. Figure 3 
shows the result of the fitting with three free parame- 
ters: An, A22, and A12. The partial density of states 
and the parameter 7 were taken from the the litera- 
ture: the two-band mapping of the four-band MgB 2 gives 
ni = 0.44 and the Fermi velocities (w 2 h )i = 3.3 and 
{v 2 ab )2 = 2.3 x 10 15 cm 2 /s 2 . 16 ' 17 This fit requires solv- 
ing two coupled nonlinear equations, Eqs. (20). We used 
Matlab with the Optimization toolbox and utilized a 
nonlinear solver using direct Nelder-Mead simplex search 
method. 20 The result is shown in Fig. 3 with the best fit 
parameters listed in the caption. We show below that 
the same set of parameters used to calculate the free en- 
ergy and the specific heat reproduces the data on C(T) 
remarkably well. 

Our numerical experimentation shows that if, in addi- 
tion to coupling constants, the partial densities of states 
and the Fermi velocities are unknown, the fitting proce- 
dure becomes unstable and good fits can be found for 
various combinations of the coupling constants. In the 
case of MgB 2 these quanities are known and our fitting 
is quite certain. 

For VsSi we do not have detailed information regard- 
ing the band structure, partial densities of states, and 
Fermi velocities on separate sheets of the Fermi surface 
of this material. Hence, we took all those as free pa- 
rameters in the fitting procedure. The conclusions thus 
are less reliable for this material than for MgB 2 : being 
mapped onto a two-band model, VsSi comes out to have 
two nearly decoupled bands with an extremely weak in- 
terband coupling (still sufficient to give a single T c ). The 
results and the best-fit parameters are given in the cap- 
tion to Fig. 4. Note that the long linear tail in p(t) as T 
approaches T c is a direct manifestation of a very small 
gap, in this case Ai, in this temperature domain. 
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with the Eilenberger expression for the energy difference: 



FIG. 3: (Color online) The data and fits of the superfluid den- 
sity for MgE>2 single crystal and corresponding temperature- 
dependent gaps (inset). The fitting parameters are: An = 
0.23, A22 = 0.08, A12 = 0.06; the partial density of states 
Tii = 0.44 and 7 = 0.56 were fixed. 




FIG. 4: (Color online) The data and fits of the superfluid den- 
sity for V3S1 single crystal and corresponding temperature- 
dependent gaps (inset). The fitting parameters are: An = 
0.1 , A22 = 0.1, A12 = 1 x 10 -5 , m = 0.47 and 7 = 0.4. 



III. FREE ENERGY AND SPECIFIC HEAT 



By fitting the data for p(t), we can extract the cou- 
pling constants X ufl along with A„(T). This allows one 
to determine all thermodynamic properties of the mate- 
rial in question, of which we consider here the specific 
heat C(T)and its the jump at T c . To this end, one starts 



F — F 

1 n ± s 

W) ' 

Near T c , one obtains: 
F — F 



(23) 



N(Q) 



7C(3)7r 2 T c 2 ^n^. 



(24) 



Following Ref . 7, one can look for solutions 8 V (t) of 
Eq. (20) near T c as an expansion: 



K = a,y 2 + b v r 3 / 2 



T=l-t. 



(25) 



We substitute this in Eq. (20) and compare terms of dif- 
ferent powers of r; the quantity — 7((3)a 2 r/2 + 
0(t 2 ). In the lowest order we obtain the system of linear 
homogeneous equations for a\i that coincides with the 
system (11). The same arguments that led to Eq. (16) 
provide 



0-2 = a-iG = a\ 



A — niAi 



(26) 



The conditions for existence of non-trivial solutions for 
b u in the next oder provide the second relation for a„. 
We omit a cumbersome algebra and give the result: 



A 2 



nin 2 rj 



1 7 C(3) A(niAn + n 2 \i 2 G 3 ) - n x n 2 r\ 
We now obtain the energy near T c , 

F n -F s = 7C(3)iV(0K 2 T 2 n»<T 2 = Br 2 

and the specific heat: 

C s -C n = ^- = UC(3)ir 2 N(0)T c £ n v al . 

The relative jump at T c is: 

AC 12 ( j. 

— n = 7W) {ni + n2G) 

A 2 — n\n 2 r\ 



(27) 



(28) 



(29) 



A(niAn + n 2 X 12 G 3 ) - rnn 2 r) 



(30) 



If all coupling constants are the same, r\ — 0, G — 1, and 
AC/C n = 12/7C(3) = 1.43, as is should be. We note 
that the sign of the interband coupling A12 has no effect 
on the jump AC since in Eq. (30) A12G is insensitive to 
this sign. 

It is easy to study numerically the jump dependence 
on the the three coupling parameters. As an example 
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FIG. 5: (Color online) Dependence of the specific heat jump 
on the mismatch between An and \22- AC/Cn at T c nor- 
malized on the BCS value of 1.43 is calculated using Eq. (30) 
with An = 0.3, A12 = 0.06 and two values of n\ and plotted 

VS A22- 
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FIG. 6: (Color online) Normalized electronic specific heat, 
C(T)/C N {T). Symbols are the data from Ref. 9. The solid 
line shows result of the 7-model, it is calculated using the 
energy (23) with parameters determined from the fit to the 
data on superfluid density. 



we show in Fig. 5 the jump dependence on the mismatch 
between An and A22 for a fixed A12 for two values of n\. 
One can see that for equal relative densities of states, the 
jump peaks at A22 = An at the one-band BCS value of 
1.43; with the fixed An and changing A22 the jump drops 
with the mismatch | A22 — An|. The peak position and 
the drop speed vary with varying bands contributions, so 
that the value of the jump per se cannot be interpreted 
as evidence for a particular order parameter. 



Now we can test our theory by employing the parame- 
ters obtained from the fit to the data on superfluid den- 
sity as shown in Fig. 3 and calculate the free energy and 
the specific heat. The result is shown in Fig. 6. The 
dashed line shows a single-gap weak-coupling BCS re- 
sults and the solid line is the result of our calculations. 
Note that this is not a fit, but a calculation with pa- 
rameters determined in independent measurement. The 
data shown by circles are taken from Ref. 9. Since this is 
NOT a fit, the agreement with the general behavior of 
C(T) and, in particular, with value of the jump at T c is 
remarkable. 



IV. T c SUPPRESSION BY NON-MAGNETIC 
IMPURITIES 

The intraband scattering does not affect T C) so that we 
focus on the effect of interband scattering with an average 
scattering time r. Since g = 1 at T c , the Eilenberger 
equations for /1.2 in two bands read: 18 

= 2Ai - 2w/i + v 2 {f 2 - fi)/r, (31) 
= 2Ai - 2w/i + v 2 (f 2 - h)/ T , (32) 



[h = 1). This system yields: 



- A 1 (^ + ni/2r) + A 2 ?i 2 /2r) 
h - -, , (33) 

_ A 2 (w + n 2 /2r) + Aim/2r) . . 

h - ; ; K 6i ) 

uicv 

where u>' = lu + l/2r. The self-consistency equation, 

A „ = ^npXvufn , (35) 

again reduces to a system of linear and homogeneous 
equations for Ai^, the determinant of which must be 
zero. Omitting the algebra, we give the result: 

P 2 n 1 n 2 r\ - P{yi\\\\ + n 2 \ 22 - nin 2 r)Q) 

+ l-Q{n l r tl +n 2 ri 2 ) = 1 (36) 

771 = nrAn + n 2 \i 2 , 772 = «2A22 + «rAi2 , (37) 

77 is defined in Eq. (12). The quantities P, Q are given by: 



y2^ = fl pa. 



2r ^ luuj' Y \9 9f 



2 ' 2tJ-H2^^ 



where t = T c /T c0 with T c0 being the critical temperature 
of the clean material given in Eq. (14). The scattering 
parameter 



1 



Po 



2irT cQ T 



(40) 
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One can easily rearrange P to the form: 



P 



1 
I 



\nt-Q. 



(41) 



Next, one solves the quadratic Eq. (36) for P and chooses 
the smaller of two roots (with the minus sign in front of 
the square root). Denoting this root as P r (X, po, t) where 
A stands for the set of all coupling constants and of partial 
densities of states, we obtain an implicit equation for t(p) 
that can be solved numerically: 



^-lnt-Q(t,p) = P r (X,p, t). 
X 

For the case rj = 0, the only root of Eq. (36) is 
1 - Q{niT)i + n 2 r) 2 ) 



(42) 



(43) 



niAn + n 2 \22 

Since in this particular case A = n\ \\\+n 2 \22 1 we obtain 
nf An + 2nin 2 Ai2 + n\X 2 2 



Ini = Q 1 



niAn + n 2 A 2 2 



(44) 



One can verify that this coincides with the suppression 
formula obtained within the model with factorizable cou- 
pling potential, see Ref. 15 or a more general work by 
Openov. 19 With the parentheses on the right-hand side 
equal to 1, this is just the Abrikosov-Gor'kov result for 
the T c suppression by a pair-braker with the scattering 
parameter po- Thus, only if n\X\\+2n\n2X\2+n\X22 < 0, 
or 



A12 < 



n\\ 2 2 



2nin 2 



(45) 



T c drops to zero at a finite r. Otherwise T c (po) is a 
decreasing function at all po. 

One can show numerically that these features of the 
T c suppression are qualitatively the same for a general 
two-band case: unlike formulas of preceding sections de- 
scribing clean materials, the sign of the interband cou- 
pling A12 does affect the T c suppression. One can verify 
that the interband scattering causes faster decrease of T c 
if the interband coupling is repulsive, A12 < 0. 

To illustrate this point we calculate suppression of T c 
with the coupling parameters An = 0.3, A22 = 0.2, and 
with three different values of A12 shown in Figure 7. 
Whereas with positive A12 the suppression is weak and 
similar to the case of materials having one anisotropic 
gap, the suppression for A12 < is much stronger. This 
finding can be checked experimentally and, in fact, the 
published data on unusually fast suppression of T c with 
carbon, aluminum, or lithium doping 21 imply that MgB2 
might have repulsive interband coupling. Otherwise, it 
is hard to reconcile the T c suppression by a factor of 
4 by 15% of C substitution. If one interprets this ef- 
fect as caused by impurities scattering, Eq. (42) with 
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FIG. 7: (Color online) Suppression of the critical temperature 
by non-magnetic impurities in a two-band superconductor. 
The intraband couplings An = 0.3, A22 = 0.2 are the same 
for all curves. The suppression is much stronger when the sign 
of A12 changes to negative (repulsive). For a strong repulsive 
A12, T c turns zero at a finite scattering parameter po- 



A12 = +0.06 provides po ~ 10 3 needed for such a sup- 
pression. This value of the scattering parameter corre- 
sponds to unrealistically short men-free path i ~ 1 Aor 
less. Changing the sign of A12 to negative (i.e. taking 
the interband coupling as repulsive), results in po ~ 0.37 
and a reasonable estimate of i ~ 400 A. The negative in- 
terband coupling would imply opposite signs of the order 
parameter on the two effective bands of MgB2, i.e. ±s 
order parameter, a proposition calling for more studies. 



V. SUMMARY 

We have presented a two-band weak-coupling 7-model 
that takes into account self-consistently all relevant cou- 
pling constants to evaluate temperature dependencies of 
the two gaps, of the superfluid density, and of the specific 
heat in clean s-wave materials. The interband coupling is 
shown to have a strong effect an these dependencies irre- 
spective of the sign of this coupling. In particular, if the 
interband coupling is negative (repulsive) it may cause 
the two order parameters to have opposite signs, i.e. the 
order parameter may have the ±s structure. In this case, 
the T c suppression by interband scattering should be very 
strong, the feature that can be utilized as a signature of 
the ±s order parameter. We speculate that a strong T c 
suppression by various dopands in MgB2 may signal such 
a possibility. All these features make the model advan- 
tageous to the empiric and not self- consistent a- model 
commonly employed to interpret the data on penetration 
depth and specific heat of two-gap materials. 
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